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Abstract We introduce a measurement scheme that 
utilizes a single ion as a local field probe. The ion is 
confined in a segmented Paul trap and shuttled around 
to reach different probing sites. By the use of a sin- 
gle atom probe, it becomes possible characterizing fields 
with spatial resolution of a few nm within an extensive 
region of millimeters. We demonstrate the scheme by 
accurately investigating the electric fields providing the 
confinement for the ion. For this we present all theoret- 
ical and practical methods necessary to generate these 
potentials. We find sub-percent agreement between mea- 
sured and calculated electric field values. 



1 Introduction 

The concept to realize and investigate quantum pro- 
cesses with single ions in Paul traps was subject to an 
impressive development over the last years. As for quan- 
tum information with single to a few ions, there has 
hardly been an experimental challenge that could with- 
stand this development, consider for instance high fi- 
delity quantum gates multi-qubit entanglement 
[3] , and error correction |4] . There has been a recent de- 
velopment in ion trap technology that was primarily ini- 
tiated by the necessity - and possibility - to scale up 
the established techniques to process larger numbers of 
ions. This idea how to scale up an ion trap computer is 
to create multiple trap potentials, each keeping only a 
small, processable number of ions [6 . Time dependent 
potentials can then be used to shuttle the ions between 
different sections. To generate such spatially separated 
potentials, the respective trap electrodes are segmented 
along the transport direction and supplied individually 
with time dependent voltages. There have also been pro- 
posals to incorporate the shuttling process into the quan- 
tum gate time by letting the ion cross laser beams or 
magnetic field arrangements [5^. Another application of 
the shuttling ability is the exact positioning of the ion in 



the mode volume of an optical cavity to perform cavity 
quantum electrodynamic experiments. Recently, multi- 
trap configurations with locally controlled, individual 
trap frequencies were proposed to generate cluster states 
in linear ion traps to perform one-way quantum compu- 
tations [TS). 

This paper presents the implementation of a method 
using a single trapped ion to investigate local electric 
fields. The method profits from the outstanding accu- 
racy of spectroscopic frequency measurements and relies 
on the ability to shuttle the ion in order to reach different 
probing sites. First, we describe the experimental system 
and the way the electric trapping and shuttling poten- 
tials are generated. After that, the measurement scheme 
using the ion as a local field probe is presented. Then, 
the quantitative experimental results are discussed and 
compared to numerical simulations of the trapping fields. 



2 Experimental field generation 

We confine a single, laser cooled ''"Ca'^-ion in a linear 
Paul trap. The radial confinement is generated by a har- 
monic radio-frequency pseudo potential resulting from a 
25 MHz drive with an amplitude of about 350 Vpp. This 
results in radial trap frequencies around 3.5 MHz. For 
the axial confinement, 32 pairs of electrodes are avail- 
able. Each of these 64 segments can be individually bi- 
ased to an voltage Vi, \Vi\ < 10 V. In the trap region the 
experiments were performed in, the segments are 250 /im 
wide and separated by 30 fim gaps (details can be found 
in [12]). The voltage data preprocessed by a computer 
and transmitted to a home-built electronics device. This 
device houses an array of 16 serial input digital-to-analog 
converters (dac, 50 MHz maximum clock frequency), 
each of which supplying four individual voltage lines, 
so that a data bus conveying 16 bit in parallel suffices 
to obtain the desired 64 outputs. The amplitude reso- 
lution of the converters is 16 bit with a measured noise 
level well below the last significant bit. The serial input 
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design of the converters keeps the amount of data being 
transmitted in parallel at a minimum. This makes the 
design extendable to a large number of channels. 

In this work, we restrict ourselves to the investigation 
of harmonic potentials. This is rectified by the fact that 
a resting, laser cooled ion only experiences the very min- 
imum of the potential. Such potentials are generated by 
applying control voltages to the trap segments. Thereby, 
each specific voltage configuration results in an axial trap 
with a given trap frequency; the position of the potential 
minimum - this is where the ion resides - can be changed 
to shuttle the ion between different locations on the trap 
axis. Doing so, two requirements have to be met by the 
time dependent voltages: (i) The amplitudes must range 
within the limits given by the dac electronics, ±10 V 
in our case, (ii) Each segment's voltage is to change as 
continuous as possible to avoid unfeasibly high frequen- 
cies. All these requirements are fulfilled by the field sim- 
ulation and voltage calculation techniques presented in 
sections [A] and [Bj respectively. For the experiments pre- 
sented here, it is advantageous to perform slow, adiabatic 
potential changes, and thus avoiding excess oscillations 
of the ion. For that reason, the update rate of digital to 
voltage converters was chosen with ~ 3.3 ms, and a par- 
allel port connection between computer and dac device 
was used. In an alternative operation modus of the trap 
supply device, the data transmission and updated rate is 
significantly sped up by using a field programmable gate 
array. Then we reach 4 fis even when updating all 64 
channels simultaneously. For the application presented 
here, it is sufficient to use a Doppler cooled probe ion. 
The probe's spatial extent can then be estimated from 
the ion temperature and the consequential extent of the 
motional wave function to be ~ 60 nm. This value could 
be reduced to ^ 5 nm by sub-Doppler cooling and a 
tighter confinement of the ion. 

3 The local field probe measurement scheme 

To test both the numerical methods and the developed 
electronic devices for accuracy, we employ a single ion 
as a probe for the electrostatic potential. We exploit the 
fact that it is possible to exactly measure the trap fre- 
quency of the confining axial potential by spectroscopic 
means. Consider, for example, a voltage configuration 
{Vi}'"'^''^ to be tested for; this configuration is meant to 
result in a trap at position Zp with frequency uj{zp). 
Then we can use time dependent potentials to shuttle 
the ion from a loading position to Zp |13| . whereby the 
final voltage configuration must be {Vi}^'^''\ The spec- 
troscopically obtained trap frequency at Zp can then be 
compared to the theoretically expected one, ujsiui{zp). 

The trap frequency is determined by measuring the 
difference frequency of the carrier and the first motional 
red sideband (rsb) excitation of the ion in a resolved side- 
band regime. For that, we excite the quadrupole tran- 
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Fig. 1 Illustration of the measurement procedure: (i) cooling 
and preparatio of the ion at starting position zq and (ii) trans- 
port of the ion to the probing position Zp. A spectroscopy 
pulse with a certain detuning excites the ion there (iii). Af- 
ter shuttling the ion back (iv), the quantum state of the ion 
is read out (v). After many repetitions for different detun- 
ings an excitation spectrum of the ion at Zp is obtained, from 
which the trap frequency (jj{zp) can be deduced. 

sition (4S'i/2,mj = -fl/2) o {30^/2, mj = +5/2) at 
729 nm. As the carrier resonance frequency does not de- 
pend on the trap frequencjj^ it is sufficient to measure 
the rsb frequency at different locations Zp. 

The measurement scheme is illustrated in figure [l] It 
consists of the following steps: 

(i) An initial voltage configuration is chosen to trap 
and cool the ion at the starting position zq. All lasers 
necessary for cooling, repumping, state preparation and 
detection are aligned to interact with the ion here. Addi- 
tionally, Zq is the position where fluorescence emitted by 
the ion can be detected by a photomultiplier tube and a 
camera. 

(ii) The ion is shuttled to the probing position Zp. 
This is accomplished by applying voltage configurations 
resulting in a series of harmonic potentials whose mini- 
mum positions lead the ion from zq to Zp. For simplicity, 



^ Nevertheless, the carrier frequency was monitored to ex- 
clude drifts in laser frequency or Zeeman splitting. 
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Fig. 2 Resonances of the first red sideband excitation on 
the S\/2 — >■ £'5/2 transition. The two peaks are measured at 
different trap positions along the trap axis. The resonance 
frequencies are given with respect to the carrier resonance 
/car, so that the trap frequency at the respective trap posi- 
tion can be read off. The full width at half maximum of the 
peaks is 8 kHz and determines the measurement accuracy. 
The different peak heights stem from slightly different field 
intensities experienced by the ion. 

these potentials were chosen such that their harmonic 
frequencies were almost constant and close to u){zp). 

(iii) Applying a spectroscopy pulse at the probing 
position Zp. Now, the voltages are exactly {Vi\'^^'\ Rest- 
ing at Zp, the ion is exposed to a spectroscopy pulse of 
fixed duration (100 ^s). The pulse is detuned from the 
(carrier) resonance by Af . This excites the ion with a 
probability P into the upper state \D^/2)- 

(iv) Shuttle the ion back to Zq, inverting step (ii). 

(v) Having arrived back at the starting position, the 
state of the ion is read out by illuminating it on the cool- 
ing transition. Whenever we detect a fluorescence level 
above a certain threshold, the ion is found in the ground 
state IS'1/2), while a low fluorescence level indicates that 
the ion has been excited to the state \Dr^/2)- 

The excitation probability for a specific detuning, 
P{Af), at the remote position Zp is obtained by averag- 
ing over many repetitions of steps (i) to (v) . By varying 
Af, a spectrum of the quadrupole excitation at the re- 
mote position is obtained without moving any lasers or 
imaging optics but the spectroscopy laser used in step 
(iii). This is a significant advantage, because that way 
even such trap positions can be investigated that cannot 
be directly observed by the imaging devices or reached 
by all lasers. 

The frequency difference between the red sideband 
and the carrier transition yields the trap frequency. Fig- 
ure[2]shows two rsb resonance peaks obtained at different 
positions within the trap volume. 

It is noticeable that each iteration cycle (i-v) includ- 
ing two ion transports results - due to the binary nature 



of the projective readout - in exactly one bit information 
about the spectrum. Therefore, thousands of transports, 
each relying on the calculated potentials, are performed 
for the determination of one frequency uj{zp). The trans- 
port, however, can be performed so fast (^ 100 /is) that 
its contribution to the overall experiment duration is sec- 
ondary; this duration is still dominated by cooling and 
detection times (~ ms). 

4 Implementation and results 

To implement the measurement scheme, we first calcu- 
lated voltage sets where z covers the whole ex- 
tent of the trap in steps of 5 /im (For arbitrary positions, 
the voltages can be interpolated). Each set results in a 
certain, wanted frequency Usimiz). That means that for 
each arbitrary position z in the trap, there can be found 
a set of voltages resulting in a potential with its min- 
imum at z and with trap frequency ujsimiz). Then, in 
order to shuttle the ion, we simply subsequently apply 
the voltage configurations for z = zq ... Zp. 

The calculated voltages are tested with high axial 
resolution, i.e. in small steps of Zp, in two far distant re- 
gions of the trap. Doing this, both small local deviations 
should be detectable and the stability over the whole 
trap structure can be tested for. To see a variation in 
uj{zp) when increasing Zp, it is of advantage that small 
variations of the trap frequency around its means value 
occur. This is a reliable way to ensure that the ion in 
fact probes the remote position. 

Figure [3] shows the expected, simulated trap frequen- 
cies and the measured ones. The data are in excellent 
agreement with the predicted frequencies. On both ends 
of the investigated trap structure, the predicted course 
o{ujsim{z) is confirmed within the spectroscopic accuracy 
of 0.6%. The mean deviation of all measured data points 
is only 0.73%. Note that the solid line shown in figure [s] 
is based solely on geometric data from a technical draw- 
ing of the trap; there is no free parameter that was used 
to match the simulations with the measurement. 

5 Discussion and outlook 

The high accordance between the simulated and the real 
trap frequencies, as can be observed in the presented ex- 
periments, indicates the reliability of at least four inde- 
pendent contributions: First, the numerical field simula- 
tion is accurate. This simulation is independent from a 
specific set of voltages, but relies on a model of the three 
dimensional trap design. 

Second, the physical realization of this geometry is 
very good, i.e. the manufacturing and assembling pro- 
cess of the micro trap is so precise that it is not lim- 
iting the field accuracy. This aspect is coming more an 
more into the experimental focus, since the miniaturiza- 
tion of the trap layouts still proceeds. For example, a 
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Fig. 3 (a) Harmonic trap frequency as a function of the po- 
sition along the trap axis. The electrode voltages were calcu- 
lated and applied such that the trap frequency shows a small 
oscillatory variation around a mean value of about 1.4 MHz 
to see the ion proceed along the axis. The solid line shows 
the frequency of the wanted harmonic potential. The data 
points are the spectroscopically measured, real trap frequen- 
cies. The error bars show the uncertainty due to the finite 
resonance linewidth. (b) Relative deviation of the measured 
from the simulated frequency \lj{z) — iUsiin{z)\ / ijJsim{z) ■ 

possible axial misalignment Ax of the two dc electrode 
wings with respect to each other would mainly result 
in an axial shift of the oscillations depicted in fig. [3^ 
by approximately Ax/2. Additionally, the amplitude of 
these oscillations decreases with increasing Ax. Indeed, 
the residuals (fig. [sja) show a small oscillatory behavior. 
This oscillation signal is shifted by at most ^15 fim with 
respect to the trap frequency oscillation. Such an effect 
could stem from a small misalignment of the two trap 
layers (see ref. |12j) or from a small offset of the theoreti- 
cal and the real, axial zero point of position. We attribute 
the observed shift to an offset error, since in this case, 
the amplitude of the oscillation is not decreased, as it 
is the case in the measured datgj^ This shift is system- 
atic and could be corrected for. Third, the calculation of 
the voltages leading to the wanted potentials is highly 
reliable. This is also qualitatively confirmed by the suc- 
cessful shuttling process itself. From our measurement 

^ The misalignment error can also be verified to be smaller 
than ~ 10 /xm from microscope pictures of the trap. 



we find that the generation of the voltages and their 
supply to the trap electrodes is indeed highly accurate. 
We also conclude that trapping fields from the applied 
voltages are not significantly perturbed by background 
fields from stray charges. 

As the generation of trapping fields in the multi- 
segment trap was found to be accurate over the whole 
considered volume, it is possible to tailor potentials for 
special purposes. For quantum computing tasks, a har- 
monic potential of constant frequency is mostly required. 
Simulations show that it is possible to obtain fixed fre- 
quency potentials with a relative frequency deviation on 
the order of 10^*. The fast transport of quantum in- 
formation in harmonic wells also requires a high degree 
of control over the trap voltages [T7]. The application 
of optimal control methods can help optimizing tailored 
time dependent potentials. Alternatively, one can feed- 
back information, gained from the measurement pre- 
sented here, into the generation process to refine the 
results iteratively jlOj. 

However, the scope of possible applications of the 
described techniques is much wider than quantum com- 
puting tasks: Ions in time dependent potentials were pro- 
posed to be used as a testbed for quantum thermody- 
namic processes [2] or for quantum simulation p31ll6| . 
The presented method is just a first proof of principle of 
using single atoms as ultra-sensitive probes. The mea- 
surement principle is not restricted to the investigation 
of the electric trapping fields. Only slight modifications 
allow for the detection of any electric or magnetic field 
with a nano-meter scaled quantum probe, as an alter- 
native to cold atom sensors [S^f^ . Then, a single atom 
field probe might be used as well for probing magnetic 
micro-structures with a relative accuracy better than 
10"'^. Also, a study of decoherence and heating effects 
and their dependence on the ion-electrode separation or 
ion location is in reach. 

Financial support from the DFG within the SFB/TRR- 
21, the German-Israel Science Foundation and by the 
European Commission within MICROTRAP, EMALI 
and SCALA is acknowledged. 

A Field calcuations 

To create a specific electric potential (/)(z) on the trap 
axis, it is necessary to find the right set of voltages 
{Vi} being applied to the trap electrodes labeled hy i = 
1, . . . ,N. The potential generated by such a set of volt- 
ages is the linear superposition of the N individual elec- 
trodes, whereby the contribution of each electrode i is 
weighted by the applied voltage Vi . After subdividing the 
axial position into a grid of M points Zj,j = 1, . . . ,M, 
we can write the overall potential at any Zj by 

N 
i=l 
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Here we introduced the electrode potential matrix Aij . It 
describes the influence of the z-th electrode to the overall 
potential at Xj . Each row i of A can be seen as a position- 
dependent function, describing the potential generated 
by electrode i (in units oi Vi), when all other electrodes 
are set to zero voltage. This quantity A is independent 
from the specific voltage and is solely given by the trap 
geometry, i.e. the shape and size of the electrode and its 
distance from Zj, for instance. 

Therewith, the potential generation can be logically 
divided into two parts: The matrix A can be calculated 
independently from voltage constraints and independent 
from the desired potential. Second, for each desired po- 
tential 0, there has to be found a set of voltages v by 
inverting the matrix equation above. Tackling the first 
problem, one recognizes that modern segmented trap ge- 
ometries can be realized in such a geometric complexity 
that conventional simulation techniques like the finite 
element method (FEM) fail. Instead, we obtained the 
potentials by solving the boundary element problem of 
the segmented trap design. Details can be found in ^iQj 

EES]. 

B Calculation of the voltages 

In the following, we address the problem how to obtain 
a set of voltages {Vi} that generates a given potential (p 
when applied to the respective electrodes. The problem 
is formally solved by matrix inversion as A~^4>. Several 
circumstances make this straight forward ansatz unfeasi- 
ble: First, there is often no exact solution to the problem, 
because is not an exactly realizable potential (note 
that in general, M ^ N is possible). In this case, an 
approximate solution has to be found. Additionally, as 
a specific electrode does not effectively contribute to the 
potential at a far distant point, this electrode's voltage 
is ill-determined. These cases have to be treated ade- 
quately by the algorithm. We solved the inversion prob- 
lem with a singular- value decomposition of the matrix A 
to identify its critical, singular values. The real N x M 
matrix A is decomposed into the product 

A = USW^, (1) 

of the unitary matrices U (N x N) and W {M x M) and 
the diagonal N x M matrix S with non-negative entries 
Sfc, fc = 1, . . . , min(Af, N). This decomposition is part of 
many standard numerical libraries and can be performed 
for any input matrix A. The inverse can then be written 
as 

A-^ = WS-^U'^. (2) 

This step is numerically trivial, because the unitary ma- 
trices are simply transposed and the entries of are 
given by 1/sfc. At this point, the advantage of the de- 
composition becomes obvious, since small values of Sk 
indicate an (almost) singular, critical value. One way to 



overcome these singular values would be to cut off their 
diverging inverse values. Instead, the Tikhonov regular- 
ization ^20^ method implies a more steady behavior as 
it makes the displacement 1/sfe — > Sk/{s\ + a^). The 
latter expression behaves like the original l/s^ for large 
values Sfc 3> a, has its maximum at Sk — a and tends 
to zero for small, critical values Sfc ^ a. From this can 
be seen that the choice of the optimization parameter a 
is a compromise between exactness and boundedness of 
the results. For a = 0, the exact solution (if existent) 
is obtained, whereas large values of a guarantee small 
inverse values and thus bounded voltage results. There- 
fore, we label the regularized quantities with index a. 
The approximate solution — (Vi, . . . , Vn) is 

v„ = WS-^U^ct>, (3) 

with being the regularized matrix with entries s^/ {s\-\- 
a2). 

Before the problem of finding an optimal a is ad- 
dressed, another constraint regarding time dependent 
voltages, i.e. series of voltage configurations, has to be 
accounted for. While moving the ion by one step, i.e. 
from a position Zp to z^, the control voltage should vary 
as less as possible. This is achieved by extending eqn. ([3]), 

= v„ + WD^W^^^, (4) 

where v represents any previous voltage set, providing 
trapping at Zp. The second term in eqn. Bl) contains a 
diagonal matrix D with entries dk = /{si. + a^). dk 
tends to zero for s/j 3> a, so that uncritical voltages are 
affected only little by the second term. For all critical 
voltages indicated by a value Sk << a, however, the 
first term in eqn. Q vanishes due to the regularization 
replacement and what remains is the contribution from 
vo, since then, d^ ~ I. Here, the choice of a determines, 
how strong the algorithm tries to generate similar volt- 
ages in a (time) series of voltage sets. 

The algorithm described above minimizes ||Av'^ — 
(/)|p-|-Q;||v^ — v|p with respect to the Euclidian norm. 
That is, the potential (p is reproduced as good as pos- 
sible under the constraint that solutions similar to the 
previous one are preferred. What remains is to find the 
proper value of a. Hereby, one has to make a tradeoff be- 
tween the boundedness of the voltages and their continu- 
ity. Under practical circumstances requiring | t^i | < Knax 
for some maximal voltage T^nax, ct can be iteratively in- 
creased to fulfill this constraint on the one hand, and 
to obtain as continuous voltage sets as possible, on the 
other hand. 
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